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ABSTRACT 

We present OH and H I Zeeman observations of the NGC 6334 complex taken with 
the Very Large Array. The OH absorption profiles associated with the complex are 
relatively narrow (Avfwhm ~ 3 km s _1 ) and single-peaked over most of the sources. 
The H I absorption profiles contain several blended velocity components. One of the 
compact continuum sources in the complex (source A) has a bipolar morphology. The 
OH absorption profiles toward this source display a gradient in velocity from the northern 
continuum lobe to the southern continuum lobe; this velocity gradient likely indicates a 
bipolar outflow of molecular gas from the central regions to the northern and southern 
lobes. Magnetic fields of the order of 200 (j,G have been detected toward three discrete 
continuum sources in the complex. Virial estimates suggest that the detected magnetic 
fields in these sources are of the same order as the critical magnetic fields required to 
support the molecular clouds associated with the sources against gravitational collapse. 

Subject headings: H II regions — ISM: clouds — ISM: individual (NGC 6334) — ISM: 
kinematics and dynamics — ISM: magnetic fields — ISM: molecules 
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INTRODUCTION 



1.1. The Zeeman Effect and Magnetic Fields 

The importance of magnetic fields in the star for- 



In this paper, we adopt the convention of RCM82 for 
radio continuum sources. 

Among the many phenomena associated with star 
formation in this complex are water masers ( Moran & 



ma ;ion process is now widely acknowledged (e.g., 
VTouschovias 1987; Shu, Adams, & Lizano 1987; Hcilcs| 



Rodriguez 198C| ), OH masers (paumc fc Mutcl 1987| ), 
methanol masers (Mentcn & Batrla 1989), bipolar 



t al. 1993| ; |McKee et al. 1993| ; |Crutcher 1999, here- 
-iftcr C9"9[). The Zeeman effect in radio-frequency 



outflows in ionized and molecular gas (Dc Pree et al 



1995; Bachiller & Cernicharo 1990; Phillips & Mam 



spe ctral lines is the only viable method for mcasur- 



paso 1991), and shocked H2 emission (Straw & Hy- 



ing the strength of the magnetic field in interstellar 



land 1989a). Several authors have discussed a picture 



mo lecular clouds and star- forming regions. If the Zee- 



ma i splitting is smaller than the linewidth, only the 
line-of-sight component of the magnetic field can be 
measured; this is always the case for non-maser lines. 

In this paper we discuss high spatial resolution Zee- 
man observations in OH and H I absorption toward 
the NGC 6334 complex. In the remainder of §1, we 
review some details of the complex; §2 contains de- 
tails of the observations and reduction of the data. In 
§3, we present the results of our observations, and §4 
contains discussions about these results. 



of sequential star formation in this complex (Cheung 
1978; [Moran fc Rodriguez 1980D in which the 



et al. 



central parts of the complex are more evolved than the 
edges. This conclusion was mainly based on the distri- 
bution of indicators of star formation such as OH and 
H2O masers and CO, far-IR and 1 mm peaks. How- 
ever, RCM82 point out that the picture is definitely 
more complex than simple sequential star formation. 
For example, source E, which lies at the edge of the 
complex, is an evolved source. 

The NGC 6334 complex has also been the subject 
of extensive recent observational studies. Kraemer et 



1.2 



The NGC 6334 Complex 



NGC 6334 is a giant molecular cloud complex and - 



al. (1997, hereafter KJPB97)| and [Kracmcr (1998a. 
hereafter K98)| observed the complex in transitions of 



star-forming region at a distance of 1.7 kpc (Neckel 



1978 ). It lies about 0°.5 above the galactic plane and 
extends ^30' parallel to it. The complex has been 
the subject of extensive studies at a variety of wave- 
lengths. Rodriguez, Canto, & Moran (1982, here 



CO, CS, and NH3 and in ionized carbon and neu- 
tral oxygen ([C II] 158 /im, [O I] 145 ^m, and [O I] 
63 ^m). They found that the molecular emission in 
NGC 6334 shows a complex structure of filaments and 
bubbles, some of which are filled with photodissoci- 



aft ^r RCM82) mapped the region at 6 cm wavelength 



wit h the VIA. They found six discrete continuum 

sources that lie along a ridge of radio emission paral- 
lel to the Galactic plane. Three of these sources had 
been known from previous single dish observations 
by |Schraml fc Mezger (1969)| . RCM82 named these 
sources A to F, in order of increasing right ascension 



ated gas (Kraemer 1998b). Boreiko & Betz (1995 



(Figures 1 & 2, this paper). pVIcBreen et al. (1979) 
mapped the region with a 40—250 /im photometer and 
found six far infrared (FIR) sources. Of these, five lie 
along the ridge of radio emission; an additional weaker 
source lies further south. The FIR sources were la- 
beled I through VI. Another continuum source, I(N), 
was seen at mil limeter and sub millimctcr wavelengths 
dCheung et al.| 1978; |Gezari| 1982); I(N) lies to the 
north of FIR I and was not detected in the obser- 



vations of 


McBreen et al. (1979) 


Dickel, Dickel, & 


Wilson (1977, hereafter DDW77) 


detected four CO 



hereafter BB95) observed the NGC 6334 complex in 
[C II] 158 ^m lines. They found that the C II radi- 
ation is bright and widespread, with a general corre- 
lation between regions of intense C II emission and 
warm dust and CO radiation. 

2. OBSERVATIONS AND DATA REDUC- 
TION 

The observations were carried out with the Very 
Large Array (VLA) of the NRAO. | The H I observa- 
tions were carried out in 1993 in the CnB and DnC 
configurations. The OH observations were carried out 
in January 1996 in the CnB configuration. Table 1 
lists the important parameters of the H I and OH ob- 
servations. For H I, both circular polarizations were 
observed simultaneously. For OH, both circular po- 



peaks within an extended region of emission. The 
nomenclature in this complex has become confused. 



3 The National Radio Astronomy Observatory (NRAO) is op- 
erated by the Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation. 
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Table 1 

Observational Parameters for VLA Observations 



Parameter 


Hi 


OH 


Date 


1993 May 27 


1996 Jan. 27 




1993 May 28 






1993 Oct. 22 




Configuration 


CnB, DnC 


CnB 


R.A. of field center (B1950) 


17 h 17*™00 s .0 


17 h 17*"OO s .O 


Decl. of field center (B1950) 


-35°45'00".0 


-35°45'00".0 


Total bandwidth (MHz) 


0.78 


0.19 


No. of channels 


128 


128 


Hanning smoothing 


Yes 


No 


Channel Spacing (km s _1 ) 


1.28 


0.27 


Approx. time on source (hr) 


12.6 


4.3 


Rest frequency (MHz) 


1420.406 


1665.402 






1667.359 


FWHM of synthesized beam 


35" x 20" 


16" x 12" 


rms noise (mjy beam -1 ) 






Line channels 


10 


8 


Continuum 


9 


7 



Note. — 1'=0.5 pc, assuming distance to NGC 6334 — 1.7 kpc 



larizations and both main lines (1665 and 1667 MHz) 
were observed simultaneously. In order to mitigate 
instrumental effects, a front-end transfer switch was 
used to reverse the sense of circular polarization re- 
ceived at each telescope every 10 minutes. Further, 
for the H I observations, all calibration sources were 
observed at frequencies displaced 1 MHz above and 
below the observing frequency for the source so that 
the Galactic H I emission would not affect the calibra- 
tion. Roberts, Crutcher, & Troland (1995) describe 
very similar observational techniques applied to SI 06. 

The editing, calibration, Fourier transformation, 
deconvolution, and processing of the H I and OH 
data were carried out using the Astronomical Im- 
age Processing System (AIPS) of the NRAO. Further 
processing of the data made use of the Multichannel 
Image Reconstruction, Image Analysis and Display 
(MIRIAD) system of the Bcrkclcy-Illinois-Maryland 
Array (BIMA). 

OH masers in the field of view can have significant 
effects upon the Zeeman analysis because they are of- 
ten highly circularly polarized. Moreover, the effects 
of strong masers can be spread over the entire field 
of view in dirty maps owing to sidelobes of the syn- 
thesized beam. Successful removal of maser effects 
is only feasible if the data are very well calibrated. 
Only in this case is the actual response in the maps 



to a point source (e.g., a maser) equal to the cal- 
culated point source response function (the "dirty" 
beam). The problem of maser contamination is espe- 
cially acute in the NGC 6334 region where there are 
two strong masers in the field. The strongest maser 
(~400 Jy) coincides with source F (Fig. 2), the next 
strongest (^50 Jy) lies within the shell-like region to 
the southwest of source A (§3.1; Fig. 2). There is also 
a weak maser at 1665 MHz (<~1 Jy) on the western 
edge of source A (Fig. 4a). Fluxes given are the sum 
of right and left circular polarizations at 1665 MHz. 
To remove maser effects from the NGC 6334 OH data 
sets, we employed a several step process. First we 
performed a self calibration on the frequency chan- 
nel with the strongest maser radiation. This channel 
is at a velocity of —8.8 km/s. Then we applied this 
self calibration solution to the entire dataset. Next we 
generated an image for each channel having maser ra- 
diation, and we specified boxes around each maser in 
the image. Within each box, we identified clean com- 
ponents that represented the maser radiation. Next 
we used the AIPS task UVSUB to compute and re- 
move from the uv dataset the effects of these masers 
by using the clean components as the input model. 
Finally, we re-generated the images for each channel 
with the maser effects removed and combined these 
images with those of the other (non-maser affected) 
frequency channels. 
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a (1950) 

Fig. 1. — Contour image of the continuum sources in the 
NGC 6334 complex at 18 cm. Contour intervals are at 
0.05, 0.1, 0.2, 0.4, 0.8, 1.6 Jy beam" 1 . The beam is 12" 
x 9"; the image has been made at a higher resolution by 
uniformly weighting the data. The rms in the image is 7 
mjy beam" 1 . Note that the two weak shell-like structures 
in the extreme south of NGC 6334, shown in Fig. 2, are 
not included in this figure. 
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Fig. 2. — Grayscale radio continuum image of the 
entire NGC 6334 complex at 21 cm. The beam is 
35" x 20", the rms in the image is 9 mjy beam -1 . 
The positions of the two strongest OH masers are 
marked by the white crosses; the northeastern OH 
maser coincides with source F. The shell-like source 
about 4' to the southwest of source A is designated 
as "shell source" in this figure. 



RESULTS 



3.1. Continuum 



Figure 1 shows the 18 cm image of the discrete 
continuum sources in the NGG 6334 complex. Fig- 



urel 9. shows the 21 cm grayscale continuum image of 



the entire complex, including the low-brightness shell- 
like sources in the extreme south that are not shown in 
Fig. 1. The 18 cm image has been made at higher res- 
olution (12" x 9") by uniformly weighting the data. 
The compact radio continuum sources in this complex 
lie in a ridge that extends from the northeast to the 
southwest; there is an underlying bed of continuum 
which extends ~5' to the northwest of the ridge. The 
compact sources are labeled A— F, as described above 



(§1.2 and Fig. 1). RCM82 and 1 


lodriguez, Canto, & 


Moran (1988, hereafter RCM88 


have discussed the 



nature and morphology of these sources. RCM82 con- 
cluded that all the compact sources, with the excep- 



tion of source B, are H II regions. Moran ct al. (1990) 
have shown that B is extragalactic. For the com- 
pact sources, our estimates of the integrated fluxes 
are in reasonable agreement (within about 10%) with 
RCM82, except for source E, where we obtain a flux 
of about 6 Jy, compared to their value of 12 Jy. 
The bipolar morphology of source A (RCM 88) is 



clearly seen in the images (Fig. 1 and 4a). In their 
higher resolution 6 cm image, RCM88 also see a dis- 
continuous shell in the core of the source. The weak 
unresolved source G351.24+0.65 to the east of source 



A was first discovered by Moran et al. (1990); De Pree 



et al. (1995) concluded that this is an optically thin 
H II region. There is also a shell-like source about 
4' to the southwest of A (Fig. 1, 2; designated as 
"shell source" in Fig. 2). Parts of it coincide with the 
position of FIR source V. [lackson fc Kraemcr (1999)] 
conclude that the relative distribution of ionized, pho- 
todissociated, and molecular gas (as seen in radio con- 
tinuum, [C II] 158 /itm, and CO 2— >1 emission respec- 
tively) toward this shell-like source conforms closely 
to an idealized model of a photodissociation region 
(PDR). In our 18 cm observations, the shell is almost 
complete, except for an opening in the south. An OH 
maser lies toward the south of this source, near the 
break in the shell (Fig. 2). 

Of the other main sources in the NGC 6334 com- 
plex, source C is extended and has a nonspherical 
appearance in the 6 cm image of RCM82. There is a 
steep decrease in the continuum intensity toward the 
south of this source (Fig. 1). In both our 18 cm and 
21 cm images, there is a low intensity source about 1' 
to the south of C. Source D is extended, amorphous, 
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T 




tended FIR source VI of [McBrccn ct al. (1979)1 Since 



-5 
v 



LSR 



5 10 
(km s" 1 ) 



Fig. 3. OH 1667 MHz optical depth profile toward (a) the 

northern lobe of source A at a=17'* 16 m 57 s .8, 5=-35°51'09"; 
(b) the continuum peak of source A at a— 17 16 rrl 57 s .8, 
S— — 35°51 42 ; and (c) the southern lobe of source A at 
a=17' l 16 m 57 s .5, <5= -35°52'27". 

and roughly spherical. Source E is also extended and 
spherical. Source F is unresolved in our observations. 
RCM82 found it to have a "nozzle-like" morphology 
with a sharp decrease in emission to the west and 
a smooth extension to the east. A very strong OH 
maser is coincident with the position of this source 
(§2; Fig. 2). 

Two weak shell-like structures appear to the south 
of the star- forming ridge in NGC 6334 (Fig. 2, 
where they are marked as "low-brightness shell-like 
sources"). These two sources are coincident with the 
bright southern parts of the optical nebula of NGC 
6334 (Kraemer & Jackson 199£). Also, part of the 
western shell-like source is coincident with the ex- 



thc faint radio shells are seen as bright objects in the 
red plate of the Palomar Sky Survey, the extinction 
toward these shells must be very low. 

3.2. OH & H I Absorption 

Optical depth profiles for the OH and H I lines were 



calculated using the procedure described in Roberts 
ct al. (1995). These OH and H I profiles all contain 
a component at about +7 km s _1 that likely arises in 
a foreground cloud. A cold H I cloud was observed at 



this velocity toward the Galactic center by ffliegel & 

TaT: 



Crutcher (1972)| and |Crutcher fc Lien (1984)| and esti 
mated to lie within 150 pc of the Sun. In all discussion 
that follows, we omit consideration of this foreground 
component and concentrate on negative velocity ab- 
sorption that is associated with the NGC 6334 com- 
plex. The H I profiles generally consist of numerous 
blended velocity components extending to velocities 
as negative as —25 km s^ 1 (e.g., Fig. 6). The only 
exception to this rule is toward source C. In this case 
the H I profiles have a single relatively narrow com- 
ponent (Avfwhm ^6 km s _1 ) centered at —6.6 km 

The OH optical depth profiles have narrower less 
blended components with center velocities in the 
range —2 to —7 km s _1 . Their center velocities and 
lincwidths are summarized in Table 2. We consider 
these components to be associated with the NGC 6334 
complex based on CO and other molecular emission 
lines (DDW77; |De Pree et al] 1995). Figure 3 shows 
the profiles in OH 1667 MHz toward the core and 
lobes of source A; the profiles toward this source ex- 
hibit a north-south velocity gradient of about 3.4 km 
s _1 pc -1 . The H I absorption toward source A is sat- 
urated; velocity structure is more difficult to quantify. 
However, the H I profiles show evidence of a north- 
south velocity gradient similar to that in OH. We do 
not detect any OH absorption toward source C; how- 
ever, detectable OH absorption is seen toward the low 
intensity continuum source about 1' to the south of 
C. Considerable velocity structure is seen in OH ab- 
sorption profiles over source D west of the continuum 
peak; we do not detect any OH absorption to the east 
of the continuum peak. Toward source E, the OH ab- 
sorption profiles are single-peaked and broad toward 
the northeast of the source, and show two peaks to- 
ward the southwest. 

The OH column density was determined using the 
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relation 




17 h 17 m 00 s 16 m 54 s 
a (1950) 




-35°43'36" 1= — i 1 l_ 

17 h 17 m 32 s 30 s 28 ! 
a (1950) 



Fig. 4. (a) Grayscale plot of N(OH)/T cx toward source A 

in the NGC 6334 complex. The contours represent the naturally 
weighted 18 cm continuum, contours are at 0.1, 0.2, 0.4, 0.8, 1.6, 
3.2 Jy beam -1 . The beam is 16" x 12". The ' + ' represents the 
position at which we have measured the B-field. The 'x' represents 
the position of the ~1 Jy 1665 MHz mascr on the western edge of 
source A. (b) Similar plot to (a), but toward source D. The contours 
are at 0.1, 0.2, 0.4, 0.8 Jy beam -1 , (c) Similar plot to (a), toward 
source E. The contours are at 0.1, 0.2, 0.4, 0.8 Jy beam -1 . 



N OH /T cx = CJ r v dv cm- 2 K- 1 (1) 

where T ex is the excitation temperature of OH and 
the constant C = 4.1063 x 10 14 and 2.2785 x 10 14 
cm" 2 (km s" 1 )- 1 for the OH 16 65 and 1667 MHz 
lines respectively flCrutchcr 1977] ). Figure 4 shows 
plots of N(OH)/T ox toward sources A, D, and E in OH 
at 1667 MHz. In principle, variations in N(OH)/T ex 
may arise due to variations in the OH column den- 
sity, or due to variations in T ex , or both. However, 
in our results and discussion, we quote only varia- 
tions in N(OH)/T ox , since the excitation tempera- 
ture cannot be measured based on absorption studies 
alone. Further discussion about possible variations in 
T cx appears in §4.2. Toward A, there is an increase 
in N(OH)/T ex just north of the continuum peak in 
the core; this ridge of enhanced OH spans the source 
from east to west. Also, we observe an increase in 
N(OH) /T cx along the eastern and western boundaries 
of the northern and southern lobes of source A. This 
effect is especially striking in the southern lobe. To- 
ward source D, there is an increase in N(OH)/T ex 
going west from the continuum peak (Fig. 4b) and 
no detectable OH absorption east of the continuum 
peak. N(OH)/T ox is maximum in the northwestern 
part of this source. Toward source E, N(OH)/T ex 
appears to increase from the center to the edge over 
most of the source (Fig. 4c). However, nothing def- 
inite can be said regarding the OH column density 
at the southeastern edge, where the information has 
been blanked in many channels to remove remnant 
effects of the strong maser in source F. 

Finally, we note the absorption profiles toward 
some of the other positions in the NGC 6334 com- 
plex. Absorption in both OH and H I is seen toward 
the shell source to the southwest of A. OH absorp- 
tion is seen toward the highest intensity patches in 
the north and west of this shell. The OH lines to- 
ward these regions have a Avfwhm ~ 3 km s _1 and 
are centered near vlsr ~ — 6 km s . Broad H I ab- 
sorption profiles with several blended velocity com- 
ponents are also seen toward the two low-brightness 
shell-like sources in the extreme south of NGC 6334 
(Fig. 2). Furthermore, OH and H I absorption is also 
seen toward some of the brighter regions of continuum 
emission to the northwest of source D. 
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Table 2 

Observed features of line profiles 



OH (1667 MHz)( a ) 



R.A. 
(B1950) 



(B1950) 



V LSR 
(km s- 1 ) 



FWHM 

(km s- 1 ) 



17 /l 16 Tl 

17 h 17 i> 
l7 h 17 7> 

17^17" 
17 h 17 v 

17 h 1? n 
17 h 17 r> 



57 s .8 
10 s . 9 
08 s . 1 
23 s . 
18 s . 7 
17 s . 
30 s . 5 
29 s . 3 



-35°5l'42 // 
-35°48'30" 
-35°49'24" 
-35°46'15" 
-35°46'39" 
-35°46'5l" 
-35°42'57" 
-35°43'l5" 



-2.5 

-5.8 
-4.7, -1.4 

-3.0 

-5.5 
-6.9. -2.2 



WA11 OH profiles also show a component at positive velocities which is not shown in this table; see §3.2 

NOTE. — H I profiles show several blended velocity components between -25 and +14 km s~ 1 toward all sources, except to" 
Dmponents; one at the positive velocity mentioned in (a) above and the other is at vj^gj^ = — 6.6 km s and ^vp^yHM ^n 



OH (1665 MHz) 




t 1 / l I=(RCP+LCP) _ 


J I I I I I 


I I I I I I I I li 
Y^-V=(RCP-LCP) 




f . f\ dl/dv scaled 
U by B, oa =148 fiG 


-10 


10 
(km s" 1 ) 




Fig. 5. — (a) Stokes I (top-histo gram) and V (bottom-histogram) 
profiles in the OH 1665 MHz line toward source A at the position 
marked by a ' + ' in Fig. 4a, CK1950 — 17 h 16 m 58 s .5, and ($1950 — 
— 35°51 / 48 // . The continuous line superposed on the V spectrum 
in the lower box is the derivative of I scaled by the derived value 
of B los = 148 fiG. (b) Same as (a), but for OH 1667 MHz, and a 
derived value of Bi os — 162 fiG. 



3.3. Magnetic Fields 

To determine magnetic field strengths using the 
Zeeman effect, we fitted a numerical derivative of the 
Stokes I spectrum to the Stokes V spectrum for each 
pixel in the absorption line cube. The technique is 
described in Roberts et al. (1993). The results of 



the fits, which give the line-of-sight component of the 
magnetic field, Bi os , toward the various sources in the 
complex are listed in Table 3. We consider the re- 
sults to be significant if the derived value of Bi os is 
greater than the 3cr level. For OH, we have imposed 
a stronger condition — the results are considered sig- 
nificant only if Bios is greater than the 3cr level for 
both 1665 and 1667 MHz lines. Field values that 
we believe to be significant are shown in bold face 
type in the table. In OH, significant detections of the 
magnetic field were made toward source A. In H I, 
significant detections of the magnetic field were made 
toward source E and source D. Figure 5 shows the 
Stokes I and V profiles toward the core of source A 
(the position marked with a '+' in Fig. 4a) in the 
1665 and 1667 MHz lines of OH, together with the 
derivative of the I profile scaled by the fitted value of 
the magnetic field. Figure 6 shows the Stokes I and 
V profiles and the scaled derivative of the I profile in 
H I toward source E. By convention, a positive value 
of Bios indicates that the field is pointing away from 
the observer. 

Toward source A, the OH data provide some ten- 
tative indications of field strength variations. Field 
strengths in the 1667 MHz line increase by a factor 
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Table 3 
Magnetic Field Values 



Source 


Position 




OH 




HI 


R.A. 


Decl. 


1665 MHz 


1667 MHz 


1420 MHz 


A 


17 h 16 m 58 s .5 


-35°51'48" 


148±20 


162 ±33 


47±15 


D 


17 h 17 m 22 s .2 


-35°46'21" 


-60±46 


-69±58 


-93±13 


E 


17 h 17 m 30M 


-35°43'09" 


-263±78 a 


-340±78 a 


169±33 












-180±29 b 


NW of D 


17 h 17 m 08 s .6 


-35°42'29" 






169±55 



a from OH data convolved with 35" beam 

b from H I data convolved with a 35" circular beam for comparison with OH results 
Note. — values in boldface type are significant detections (see §3.3) 



a i 



. ' H I (1420 

~, i , , , 


MHz) 

-ijjAl^RCP+LCP) 




ft 


l-i/^VHRCP-LCP) 




, C\ dl/dv scaled by 
V B los =-169 



-50 50 

t lsr (km a" 1 ) 



Fig. 6. Stokes I (top-histo gram) and V (bottom-histogram) 

profiles in the H I line toward source E at the position 01950 — 
17 h 17*™30M, and ^1950 = -35°43'09". The continuous line su- 
perposed on the V spectrum in the lower box is the derivative of I 
scaled by the derived value of B los — —169 fiG. 

of two or more to the north of the continuum peak, 
in the direction of the east-west ridge of enhanced 
OH optical depth (Fig. 4a). However, weak maser 
contamination in this ridge at 1665 MHz prevented 
us from confirming this result in the 1665 MHz line. 
This maser coincides in velocity with one side of the 
absorption line, and its effect in the Stokes V profile 
cannot be removed completely. We can, of course, ex- 
clude the channels contaminated by the maser from 
the fit, and compare the values of the magnetic field 



in 1665 and 1667 MHz. However, higher resolution 
OH data will be needed to investigate such variations 
anyway since the existing data provide only three to 
four independent resolution elements across source A. 
Therefore, we defer examination of the spatial varia- 
tion of the magnetic field in source A to a later higher 
resolution study. 

Toward source E, where we have a significant de- 
tection in H I, we also find a possible detection in 
the 1665 and 1667 MHz absorption lines of OH. Our 
maser removal process has diminished the effect of the 
mascrs to a large extent. However, it has proved im- 
possible to completely remove the effect of the strong 
maser in source F in the Stokes V profile in a region so 
close to the maser. Therefore, in fitting for the mag- 
netic field in source E, we have fitted over only those 
channels in the absorption line that are not affected 
by the maser. To improve our noise estimates, we 
have convolved our OH data with a 35" beam. For 
comparison, we also convolved the H I data with a 
circular 35" beam. Magnetic field values in the 1665 
and 1667 MHz OH data are consistent, and the H I 
data reveal a significant field of similar order and the 
same sign. Finally, there is also a possible detection 
in H I toward an area of continuum to the northwest 
of source D. This fit is barely at the 3cr level, so we 
state it as an upper limit. If it is real, it represents 
a reversal in the sign of the field in going from the 
region near D and E to the region northwest of D. 
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DISCUSSION 



4.1. Individual Sources in the NGC 6334 Com- 



plex 



4.. 1.1. Source A 

The morphology of continuum emission from this 
source, including its bipolar structure, has been de- 
scribed in §1.1. The peak of the FIR source IV lies 
near the peak of the radio source. Also associated 



with source A are a Herbig-Haro-like object (Gyul 
audaghian, Glushkov, & Dcnisyuk 1978; Bohigae 



_ ^„~~— , — -~ , ^ . ^ p |— 

1992|), H 2 Q mas crs ( |Rodrigucz ct al. 1980| ; |Moran fc 
Rodriguez 1980|), and a CO peak (DDW77; |Phillips 



ic Vries, fc de Graauw 1986, hereafter PVG86| ; K98) 
RCM88 suggested that the bipolar morphology of this 
region results from the confinement of the H II region 
by a flattened structure of gas and dust. This east- 
west confinement allows the ionized gas to escape only 
to the north and south. Indeed, a 40" resolution map 



of 100 /im optical depth from Harvey & Gatley (1983) 



shows an east-west elongation, which suggests a flat- 
tened structure of dust with an east-west extent of 
at least 100". Their map center coincides with the 
center of the ridge of enhanced N(OH)/T cx in our 
data (§3.2 and Fig. 4a). The east- west ridge of en- 
hancement in our OH optical depth may be related to 
this disk. KJPB97 observed various transitions of CO 
and CS toward source A (HPBW of CO J=1^0 ob- 
servations = 53"). Their data show that source A is 
surrounded by a 4'. 4 x 1'.8 structure, indicating that 
the molecular cloud is also flattened on much larger 
scales. 

The kinematical data for source A and its environs 
also offer insights into the physical nature of this re- 
gion. On the larger scale of 4' (2 pc), KJPB97 find 
an east-west velocity gradient of 2.4 km s^ 1 pc -1 in 
the gas traced by CO and CS. They conclude that 
a rotating disk of molecular gas surrounds source A. 
Our OH absorption data also reveal an east-west ve- 
locity gradient of the same sign. Over a 35" (0.3 pc) 
scale size, this gradient amounts to 4 km s _1 pc -1 



Do Pree et alj (1995) have observed ionized and neu- 
tral gas (H76a and H92a recombination lines, and 
H2CO absorption) toward source A with a resolution 
of about 5". They report some indications of east- 
west velocity gradients in the opposite sense of those 
cited above. KJPB97 have noted the disagreement 
between the sense o f rotation as co nveyed by their 
molecular data, and Dc Pree et alj (1995)'s recom- 



bination line and H2CO absorption data. However, 
our OH data agree in sign with KJPB97's data over 
a larger scale size. The velocity gradients seen by Da 



Pree et al. (1995) over smaller size scales (about 15") 



may, therefore, reflect local velocity inhomogeneities 
in the core of the source, rather than any bulk rota- 
tion in the east- west direction. 

Evidence for north-south motions in the bipolar 
lobes comes from two sources. De Pree et al. (1995) 



detect a north-south velocity gradient in H92a over 
a size scale of 100", amounting to 19 km s _1 pc^ 1 . 
They conclude that this indicates a bipolar outflow 
of ionized gas into the northern and southern lobes 
from the central source. They argue that this out- 
flow lies nearly in the plane of the sky, based on the 
extended bipolar morphology and small velocity gra- 
dient. Our OH absorption data (§3.2 and Fig. 3) 
reveal a north-south gradient of 3.4 km s _1 pc -1 over 
a size scale of 80". This suggests that the molecular 
gas traced by OH absorption is also part of an outflow 
from the central regions to the northern and southern 
lobes. This north-south velocity gradient in OH has 
the same sign as the north-south gradient in H92a 

(1995). This correspon- 



dctcctcd by Dc Pree et al 



dence in sign suggests that the OH gas may lie close 
to the neutral-ionized interface and be entrained by 
the outflow of the ionized gas. 

The [C II] 158 [im profile of BB95 (HPBW = 43") 
toward FIR IV shows two peaks near —1 km s _1 and 
—5 km s _1 . BB95 note, however, that their spectra 
cannot be fitted by two Gaussian emission compo- 
nents. They mention a feature toward FIR IV with 
a velocity near —3 km s" 1 which is apparently due 
to self-absorption. We convolved our OH data with a 
43" beam for comparison with their data. Toward the 
position of FIR IV in source A, our OH profiles show 
an absorption feature near —3.7 km s~ 1 . Therefore, 
our data support the idea that the minimum between 
— 1 and —5 km s _1 in the [C II] 158 /im profile is, 
indeed, self-absorption. The CO profiles of DDW77 
and PVG86 also show a peak near —3 km s _1 toward 
FIR IV, and a Avfwhm ~ 7 km s _1 . 

A significant Bi os (~ 150 /iG) was detected in OH 
1665 and 1667 MHz absorption toward source A. We 
note that the magnetic field obtained from the formal 
fit in H I is of the same order although less than the 
magnetic field detected in OH. Due to the blending in 
velocity in the H I gas, it is almost certain that optical 
depth effects will suppress any Zeeman signature in 
the V profile ( Schwarz et al. 1986 ) , which may explain 
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why we do not have a more significant detection of 
Bios in H I toward source A. 

4-1.2. Source D 

The extended, featureless source D (Fig. 4b) coin- 



cides with the FIR peak II of McBreen et al. (1979) 



and an H2O maser (Moran & Rodriguez 1980). Based 
on their 6 cm radio continuum measurements, RCM82 
found that a ZAMS 06.5 star would be required to 



ionize this H II region. Straw, Hyland, & McGregor 



(1989) found that IRS 24, which lies close to the cen- 
ter of the continuum source and has a dereddened K 
magnitude corresponding to a ZAMS 06 star, is the 
dominant source of excitation in this source. Source 
D is located at the edge of a hole in the molecular 
gas emission, as seen in K98's CO data. Their mod- 
els suggest that the molecular hydrogen column den- 
sity toward source D is t he lowest of any source in 
the NGC 6334 complex, fctraw fc Hyland (1989b) 



used near-infrared star counts to find the extinction 
through NGC 6334; their estimates also suggest that 
the extinction toward source D is less than that to- 
ward any other sources. This hole in the molecular 
gas emission is clearly reflected in our high resolution 
image of source D (Fig. 4b), where no OH absorption 
is detected east of the continuum peak. 

BB95 model their C II profile toward FIR I (source 
F) as an emission component centered at vlsr = —5.2 
km s _1 and an absorption component near vlsr = 
— 1.6 km s — 1 , and note that the source of the absorb- 
ing gas at this velocity is unclear. This absorption 
feature is also found near FIR II (toward source D). 
The CO profiles nearest to FIR II (of DDW77 with a 
70" beam and of PVG86 with a 1'.7 beam) are broad 
(Avfwhm ~ 7 km s _1 ) but appear to be Gaussian 
and single-peaked near —3 km s _1 . Our OH absorp- 
tion profiles in this region show considerable veloc- 
ity structure (§3.2). However, when convolved with 
a 43" beam, our OH spectra show two components 
near vlsr ~ —5.0 and —1.6 km s _1 . This result sug- 
gests that both C II velocity components identified by 
BB95 lie on the near side of the H II region source D. 

4.1.3. Source E 

Source E is the northernmost compact radio con- 
tinuum source in NGC 6334. RCM82 found that a 
ZAMS 07.5 star would be required to p roduce this 
H II region. Tapia, Persi, fc Roth (1996) detected at 
least 12 faint K-band sources toward source E; these 



sources were not detected in the J- or H-bands. They 
suggested that these are a cluster of BO-BO. 5 ZAMS 
stars which collectively ionize the H II region. 

The C II profile of BB95 closest to source E is 
the one 1' northwest of FIR I. The profile displays a 
peak near ~ —8 km s , and a lower level wing ex- 
tends up to ~ +2 km s _1 . Interestingly enough, the 
OH absorption profile toward the north of source E, 
when convolved with a 43" beam for comparison with 
BB95's data, is almost identical to this C II profile in 
shape and extent, except that the peak of the absorp- 
tion is at —6.6 km s . The H I absorption is present 
at all these velocities. 

4.2. Virial Estimates 

A principal goal of Zeeman effect measurements 
is to estimate the importance of the magnetic field 
to the dynamics and evolution of star forming regions 
like NGC 6334. Such estimates also require that other 
physical parameters of the regions be known, such as 
internal velocity dispersion, proton column density, 
radius, and total mass. C99 has summarized these 
concepts. Here, we apply such a simple analysis of 
magnetic effects to NGC 6334 sources A, E, and D, 
especially source A for which physical quantities are 
best known. 

First, we must estimate other relevant physical pa- 
rameters. Column densities can be estimated from 
molecular line data such as that reported here, sub- 
ject to several assumptions. For source A, we estimate 
column densities from our own OH data. We use an 
average value of N(OH)/T ox equal to 6 x 10 14 cm~ 2 
K _1 to calculate N(OH) toward this source. The ex- 
citation temperature cannot be measured based on 
absorption studies alone. We adopt T cx = 40 K. K98 
fou nd Tk — 50 K toward source A based on CO stud- 
ies. |Forster et al. (1987)| found T k - 40 K from their 
ammonia studies. It is unlikely that T cx remains con- 
stant over this source, and it may vary down to the 
dark cloud value of 10 K in the densest parts of the 
source. However, the results will scale with the tem- 
perature. The c onversion ratio Nh?/N (OH)=2 x 10 6 



was taken from Roberts et al. (1995) for S106. This 



ratio is an order of magnitude higher than the ra- 



1979 



tio for dark clouds determined by Crutcher 
However, we find N(OH)/A v ~ 5 x 10 14 , where we 
have used A v = 50 mag., based on PVG86's estimate 
of A v ~ 40 mag. toward a position about 1' to the 
southeast of A. Using the standard conversion ratio 
of N H2 /A V ~ 10 21 , we obtain N H2 /N(OH) ~ 2 x 10 6 
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Table 4 
Source A Parameters 



Parameter Value 

Radius (r) 0.8 pc 

Mass (M) 2200 M 

N H2 4.8 x 10 22 cm- 2 

n p 3 x 10 4 cm" 3 

T k 40 K 

Av 2.6 km s' 1 

Bi os 150 fiG 



in agreement with [Roberts ct al. (1995) 



lies our use of the Roberts et al. (1995) 



which justi- 
value rather 

than the dark cloud value. Using this conversion ra- 
tio, we obtain a value of Nh 2 = 4.8 x 10 22 cm~ 2 . 
Note that this value compares well with the value of 
N H2 = 1.3 x 10 23 cm~ 2 obtained by KJPB97 based 
on their CO excitation models, and with K98's value 
of 8.9 x 10 22 cm~ 2 based on CO integrated inten- 
sity. Next, the most problematic parameter is the 
radius of the molecular cloud. K98 report a range of 
values for the radius based on their CO and CS emis- 
sion studies. Using the average value of Nh 2 =4.8 
x 10 22 cm' 2 and different values for the radius, we 
estimated the mass and proton density toward source 
A. The best value of the radius was then taken to be 
that for which the mass of the molecular cloud most 
closely matched the mass obtained by KJPB97 from 
their CO data. Further, this adopted value of the 
radius (0.8 pc) matches well with the radius of the 
molecular cloud in 13 CO 2->l (0.8 pc) and CS 3^2 
(0.7 pc) from KJPB97. Results of our analysis for 
source A are given in Table 4. Note that our value 
for the proton density is ~ 2 higher than their value. 
For source E, where some of the above parameters 
are not available, we adopt a value of Ah 2 =6.3 x 
10 22 cm~ 2 from K98, as derived from CO excitation 
model calculations. Similarly, for source D, we use 
their value of Nn 2 = 2.2 x 10 22 cm" 2 . 

The most straightforward estimate of the impor- 
tance of the magnetic field comes from the relation 



B s ,crit = 5 x 10~ 21 N p fiG 



(2) 



where N p is the average proton column density of the 



cloud, and Bs, C rit is the average static magnetic field 
in the cloud that would completely support it against 
self gravity. Note that the magnetic field in a cloud 
can be described in terms of a static component P>s, 
and a time-dependent component Bw- It is likely that 
the Zeeman effect primarily samples the static compo- 
nent (see Brogan et al. 1999| , and references therein). 
If the actual field in the cloud is comparable to Bg. cr it, 
then the field can be judged dynamically important 
to the region. For source A, we find Bs, cr it = 500 /xG, 
compared to an estimated total field strength of 300 
fiG. Following C99, we have used total magnetic field 
strength B equal to 2 times Bi os , and |B| 2 equal to 
3 times |Bi os | 2 . For source E, we find Bs, C rit = 600 
fiG, compared to an estimated total field strength of 
400 fiG. Again, for source D, we find Bs.crit = 220 
fxG, compared to an estimated total field strength of 
190 ^G. Therefore, subject to the obvious uncertain- 
ties in such estimates, we judge the magnetic field in 
each region to be comparable to if slightly less than 
the critical field. In all three regions, the magnetic 
field should be dynamically significant, providing an 
important source of support against self gravity. 

Other parameters of magnetic significance for source 
A can be inferred from the data of Table 4; they are 
shown in Table 5. All symbols used in the discus- 
sion that follows are described in Table 5. We find 
er/c =2.9, indicating that the motions in the cloud 
are supersonic. This is a common result in molecular 
clouds; C99 has found that motions are supersonic 
by about a factor of 5. However, the ratio ct/va in- 
dicates that they are sub-Alfvenic, suggesting that 
these motions may be the result of Alfven waves in 
the cloud. The parameter j3 v « 0.03, close to the 
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Table 5 

Source A : Derived Values and Virial Estimates 



Parameter Value 

cr/c 2.9 

Ct/va 0.3 

/3 p 0.03 

[M/$ B ] ofcs/c „ t 1.6 

W 3.7 x 10 47 ergs 

T/W 0.22 

M/W 0.27 



Note. — a is the line velocity dispersion, c is the speed of sound, va is the Alfven velocity, M/$b is the mass to 
magnetic flux ratio, W is the virial gravitational energy, T is the virial kinetic energy, and M. is the virial magnetic 
energy; all terms are defined in C99. 



average value of 0.04 determined by C99, suggesting 
that the magnetic fields are important, and that the 
magnetic pressure dominates over the thermal pres- 
sure. The observed mass-to-magnetic flux ratio with 
respect to the critical value is 1.6, indicating that the 
cloud is magnetically supercritical. Also, T w Ai. 
Thus, the kinetic and magnetic energy densities are 
in approximate equilibrium, which would be expected 
if the magnitude of the fluctuating part of the mag- 
netic field was approximately equal to the magnitude 
of the static field (C99). Finally, we find that 2T + 
A4 i=a 0.7W. Hence, since the external pressure term, 
which is not included here, will act in the same sense 
as the gravitational term, the cloud is in approximate 
virial equilibrium. 

5. CONCLUSIONS 

1). We have observed OH and H I in absorption 
toward NGC 6334 with the VLA. The OH absorption 
profiles are relatively narrow (Avfwhm ~ 3 km s _1 ) 
and single-peaked toward most of the sources. To- 
ward source A, the OH profiles display an east-west 
velocity gradient; this gradient has the same sign as 
the velocity gradient detected in molecular emission 
lines by KJPB97. The OH profiles toward source A 
also display a north-south velocity gradient which is in 
agreement with a similar gradient discovered in H92a 
recombination lines, and which likely indicates a bipo- 



lar outflow of material into the lobes seen in the con- 
tinuum source. No OH absorption is detected toward 
the continuum peak of source C itself. Toward source 
D, the OH profiles show considerable velocity struc- 
ture west of the continuum peak; no OH absorption is 
detected east of the continuum peak. The H I profiles 
toward NGC 6334 contain several blended velocity 
components, except toward source C, where the H I 
profile at negative velocities is almost single-peaked. 

2) . We have calculated N(OH)/T cx toward the 
continuum sources. Toward source A, we observe an 
increase in N(OH)/T ex in an east-west ridge north of 
the continuum peak; there is also an increase toward 
the eastern and western boundaries of the northern 
and southern lobes. 

3) . Magnetic fields of the order of 200 /zG have 
been detected toward some of the sources in this com- 
plex. In OH, we have significant detections of mag- 
netic fields toward source A, and possibly a significant 
detection toward source E. In H I, we have significant 
detections toward sources E and D. There may also 
be a detection toward an area of continuum to the 
northwest of source D. Also, Bi os obtained from the 
formal fit for source A in H I is of the same order as 
Bios detected in OH absorption toward source A; sim- 
ilarly, Bios obtained from the formal fit for source D 
in OH matches Bi os detected in H I absorption toward 
this source. 
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4) . Wc have three or four independent beams 
across source A; measurements at different positions 
in 1667 MHz suggest that the magnetic field increases 
going north of the continuum peak. The signs of the 
detected magnetic fields are opposite in going from 
sources E and D (where we have a detection in H I 
absorption) to source A (where we have a detection 
in OH). If the field detected in H I to the northwest 
of D is real, there is also a sign reversal in the field in 
going from D and E to the position northwest of D. 

5) . We have used various observed and derived 
parameters to study the implications of the detected 
magnetic fields toward sources A, E, and D. In all 
cases, it appears that the detected fields are of the 
order of the critical field needed to support the molec- 
ular cloud associated with that source against grav- 
itational collapse. We have also compared various 
derived parameters for source A with the published 
values of those parameters determined from an en- 
semble of clouds by C99. 
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